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We cast the four-dimensional field equations of the Nonsymmetric Gravitational Theory (ngt) 
into a form appropriate for numerical study. In doing so, we have restricted ourselves to spherically 
symmetric spacetimes, and we have kept only the Wyman sector of the theory. We investigate the 
well-posedness of the initial-value problem of NGT for a particular data set consisting of a pulse 
, in the antisymmetric field on an asymptotically fiat space background. We include some analytic 

0^ ' results on the solvability of the initial- value problem which allow us to place limits on the regions 

of the parameter space where the initial-value problem is solvable. These results are confirmed by 
. numerically solving the constraints. 

o ■ 

; I. INTRODUCTION 

^ ■ Of late, there has been an increasing interest in numerical relativity and the application of numerical techniques to 

' the study of dynamical gravitational objects. The reason for this is fairly clear: numerical relativity makes it possible 

. to study physical situations that have defied analytic analysis. 

' The Nonsymmetric Gravitational Theory (ngt) (see , as well as j7| for the most recent discussion) 

, is an alternate theory of gravity that evolved from a re-interpretation of the Einstein Unified Field Theory (uft) 

■ (see P,p|, p^ , pl] ,p^ 13 1^,|1^), itself an extension of Einstein's General Relativity (gr). The theory proposes to do away 
\ with the assumption that the metric tensor describing spacetime is a symmetric, second-rank tensor, replacing it with 

■ a second-rank, nonsymmetric tensor. This tensor is typically written as the metric tensor in GR: However, since 



X 



a metric tensor is by definition a symmetric quantity, g^i/ cannot represent the metric of the ngt spacetime; instead. 



a metric tensor is by aennition a symmetric quantity, g^i/ cannot represent tnc metric oi tnc NC 
I . we refer to it as the fundamental tensor (see the discussion in ||l^,|ri],^ and in particular fl^ ) 



Although on the surface the above-mentioned substitution appears innocuous, in reality it induces fundamental 
changes in the underlying structure of the theory, both physical and mathematical. Indeed, because of the lack of 
a definite tensorial symmetry in the various quantities of the theory, the field equations of ngt are somewhat more 
complicated than the corresponding field equations of GR, as is the dynamical structure. For instance, it can be shown 
that a spherically symmetric gravitational field in GR is diffeomorphic to a static, spherically symmetric solution, so 
d ' that all spherically symmetric stellar objects have a static exterior. This is known as Birkhoff's theorem (see for 
instance |jl^, p. 372). In contrast, it was demonstrated in that there is no equivalent to Birkhoff's theorem in 
NGT. Hence, if we wish to study a stellar object even so simple as a spherically symmetric star in NGT, we must 
consider the full, time-dependent theory. 

The increased mathematical complication of the field equations of NGT versus those of GR has meant that NGT has 
been plagued since its inception by a distinct lack of exact, or in some cases even approximate, solutions. In particular, 
in addition to the (trivial) GR solutions, a single static solution to the spherically symmetric field equations is known to 
exist for massless NGT: the so-called Wyman solution (sec |p^ , p^ , ^ , ^ ) . The so-called Papapetrou solution (see [p2[ , 
in particular SI) has been formally eliminated from the theory (see M). The addition of terms to the action by 
Moffat does not simplify matters: exact solutions for the massive theory are even more elusive than for the 

massless theory. It therefore appears that NGT would be an ideal candidate for numerical investigations. 

In this work, we begin the task of defining a formalism that will allow the application of the techniques developed 
in the field of numerical relativity to NGT; the work was begun in |2^,^,^, where a Hamiltonian formulation of 
NGT was derived from the action principle of the theory. It was found in these references that the antisymmetric 
degrees of freedom separate into two sectors, namely the Wyman sector and the Papapetrou sector. One of these 
sectors, specifically the Papapetrou sector, is plagued with linearization instabilities; this is, in fact, what motivated 
the work in 0, which eliminates these degrees of freedom. However, since the two sectors of the theory are decoupled 
in evolution, by choosing Wyman sector initial data, our results are equally applicable to most UFT-like actions. 
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The next section gives a brief overview of the coordinate choices and shcing conditions that allow us to rewrite 



the set of first-order field equations in a form appropriate for numerical study. In §111 we discuss the initial- value 



problem, as well as some difficulties that arise when dealing with a constrained set of field equations, in particular 
how those problems apply to ngt. In we discuss the numerical techniques relevant to the solution of the initial- 
value problem of NGT. We also present one such generated solution and discuss its properties. The first appendix 
contains the spherically symmetric reduction of the field quantities, while the second appendix gives the conformal 



transformation rules of the field equations. The last appendix is a brief demonstration of two results used in §111 



II. THE SPHERICALLY SYMMETRIC NGT FIELD EQUATIONS 

Early versions of not consisted of a re-interpretation of UFT, where the antisymmetric components of the 

fundamental tensor, originally introduced to represent the components of the electromagnetic field, were considered 
gravitational degrees of freedom. The work in [^,|7j necessitated the addition of two terms to the action. One of 
these, the so-called Bonnor term (see endowed the antisymmetric fields with a mass and introduced a parameter 
/i into the theory. This parameter is referred to interchangeably as the NGT mass parameter or the inverse of the 
NGT range parameter. The findings of (25| lead to a further modification of the theory, described in [Q. This left 
/i as a free parameter of the theory, with a potential phenomenological interpretation in applications such as galaxy 
dynamics (see p9|). 

Rather than re-iterating the full set of four dimensional NGT field equations here, we will simply refer the interested 
reader to the literature (sec, for instance, [P|J|,^,||,|[0] , [0, and among others). Instead, we offer in Appendices A 
and B an overview of the spherically symmetric, Wyman sector field equations written in first-order form; these will 
be sufficient for our purposes. These appendices serve also to define the notation we use in this work. The goal of the 
current section is to make the coordinate choices and to choose the slicing conditions that will reduce the NGT field 
equations to a form amenable to numerical study. 

Upon performing a conformal transformation of the field variables: 7'^'' cf)^'^^"-^ and Kab — > (f>~'^Kab, where 7"'' 
is the inverse of the hypersurface metric and K^b is its extrinsic curvature (for details, see Appendix B), we find 
that, just as in GR, the Hamiltonian and momentum constraints remain inter-twined. In order to make any headway 
in solving these constraints, we must find a way of decoupling them. In GR, this is accomplished, for example, 
by using maximal slicing: Tr[i4r] = K'^^ = l^'^K^^, = 0, which fully decouples the Hamiltonian and momentum 
constraints. Preserving this condition in time, ^tiTrl-fi^]] = 0, yields a second-order elliptic equation for the lapse 
function: 7°''VaVb[iV] - TVi?^^) ^ q (ggg^ ^j. instance, (3.4) in §IILB of fs^). A similar approach can be taken in 
NGT; however, the resultiiitf "lapse equation" is not nearly as simple as its GR counterpart. Instead, it turns out to 
be more profitable to taketl K^2 = 0. In this case, the symmetric-sector components of the extrinsic curvature drop 
out completely from the Hamiltonian constraint. Although both the conformal factor and the extrinsic curvature still 
appear in the momentum constraint, it is only the Kn component of the latter that appears. Therefore, if a solution 
to the Hamiltonian constraint can be found, this can be inserted into the momentum constraint, which can then be 
solved for the remaining component of the extrinsic curvature. 

In spherical symmetry, there is only one coordinate choice to be made: a radial coordinate r is chosen such that 
[(^22-j2 _|_ j-^[23]-j2ji/2 _ j^-^ (^j.'j ^ where R{r) is some currently unspecified function; we will choose R{r) = r. This 
choice can obviously be satisfied by taking 7^^ = R~^(r)cosip and 7'^^! = i?^^(r) sin^, or tan^/j = 7'^'^l/7^^. The 
components of the surface Ricci tensor are then (see Appendix A) 

-2d^[\nR] - 2{dr[lnR]f - ^{dr[^]f - 9^ [In 7" ] 9^ [In i?], 

1 - i?27"(cos V'<9^[lni?] -^2cos'iP{dr[lnR]f ~ i sin^pdf.[^j] ~ sin V^^^ [In i?] 9^ [V'] ) 

- ii?25r[7"](cosV5r[lni?] + i sin V'9r[V']), 

-i?^7"(sinV'5^[lni?] - 2 sin 1/; (9,. [In i?])^ + i cos ^'^^[V^] -t- cos7/'9r[lni?]a,.[^]) 

- ii?2a,.[7"](sin7/)9,,[lni?] - \cos.i\}drW\), 

from which the Ricci scalar is found to be 



and 



pNS(3) 



R. 



NS(3) 



R 



NS(3) 
[23] 



iProm Appendix A, K^2 = 7^^-ft^22 + 7'^^'i[23l and fa = 'y^^^^K22 - 7^^i[23]- 
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^NS(3) ^ _27"(2a2[lnii'] + 3(5,[lni?])2) + 2i?-2cosV' - i7"(9.[V^])2 - 29,[7"]a,[ln i?]. (1) 
We also define 

R\ = 722<'^'^ +7'''li?p3p^ = i?-' cos 7/^ - 7"9,'[lni?] - 27" (5, [In i?])^ - i9,[7"]a,,[lni?], 

and 

As mentioned above, if we clioose tlie slicing condition K^2 = On the Haniiltonian and momentum constraints 
are effectively decouple d, a llowing them to be solved independently .□ In fact, using the definition of the derivative 
operator A^[] given in ( pl[ ), the Hamiltonian constraint is written 

H = = [7"]-i/20[8A^[0] - 0i?NS(3) ^ 2i?2^-7[(/3)2 _ K\{K\ + 2K\)] + ^t^^R^cf^' sin^ V] (2a) 
= [7"]-i/2^[8i?2^"52[(/,] + 16i?7"a,[0] +4i?2a,[0]9,[7"] + ^R^^^^^drW)^ 

+ 2R(l)dr[-/^^] ~ 20COS + 2(/)7" + 2R^ {f :if (j)-"^ + I fJ-"^ R"^ (jj^ sin^ V], (2b) 

while the momentum constraint becomes Tii = — 2i?^[7^^]^^/^(7^^i4rii9,,[ln(i?^(/)^)] — j^39r['0]) — 0. The latter one 
of these is solved for the sole remaining symmetric-sector component of the extrinsic curvature in a straightforward 
manner: 

" 7iia,[ln(i?2^4)]- y-') 

That the momentum constraint can be solved for Kn in this fashion allows us to implement a "semi-constrained" 
evolution. What is more, by setting (f> = 1 in (pb|), we obtain a first-order, ordinary differential equation for 7^^, whose 
solution on every time-slice would generate the evolution of that field variable. This would allow us to implement a 
fully constrained evolution, eliminating both the Hamiltonian and momentum constraints, as well as the symmetric- 
sector evolution equations, hence by-passing the problems inherent in evolving a set of discretized field equations 
(see, for instance, the discussion of the Cauchy problem in as well as the issues raised in and on the 
maintenance of constraints during numerical evolution). However, at this time there does not seem to be any obvious 
advantage in doing this, and it would only further cloud the issue of the solvability of the Hamiltonian constraint 



(see ^HI below). In this work, we have chosen to implement the semi-constrained evolution mentioned above. 

Naturally, the slicing condition = and the coordinate choice [(7^^)^ -|- ('-j^^^^y]^/^ = R~^(r) mentioned above 
must be enforced: dt[K'^2] — and dt[R] = 0. The first of these leads to an equation for the lapse function: 

-[7"9,'[02] -f i?-«<^a.[0]a,[7"i?^] + iA^'^'sin^^]. (4) 



TV 7"9^[i?02 



This first-order, ordinary differential equation is solved on every time-slice for the lapse function. Meanwhile, pre- 
serving the coordinate condition in time leads to an equation for the only non-trivial component of the shift function: 

dr[N^] - 2N^dr[\n{R(f>^)] - 2N(p-^K^2 = 0. (5) 

When K'^2 = 0, the solution of this is = AR'^cp'^, where A{t) is a constant of integration; here, we take A{t) = 0, 
so that — > as r — !■ +00, maintaining asymptotically fiat spatial slices in evolution. 

Once the slicing conditions, coordinate choices, and constraints are taken care of, the evolution equations are: 

947"] = -2TVr'[7"]'^ii (6a) 
dt[ij] = 2N^-^j\ (6b) 
dt[Kii] = 4iV[092[0] _ (a,[0])2] + 2cbdrmNdr[lnj'^] + 2iV9„[lni?] - dr[N]) + cj}^d^[N] 

+ ^cb^dr[N]dr[\nj^^] - N(f>^R^^^^'^ + N^-'^j^^iKn? + i/i'^^'b"]"' sin^ (6c) 
dtifs] = 7V07"<9,[</)]c),[V'] + i(/.27"c),[iV]a,[i^] - Ncj^^R^ - Ncl>-''-f^^ Knfs - ^fi^Ncj)'' sin./, cos V'. (6d) 
Note that we have included here the evolution equation for Ku, for the sake of completeness. 



^Strictly speaking, the Hamiltonian constraint must be solved first, and its solution must be used to solve the momentum 
constraint. 
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III. THE INITIAL- VALUE PROBLEM IN NGT FOR SPHERICALLY SYMMETRIC SYSTEMS 



In the words of York and Piran, "The initial value problem of a physical theory is ascertaining what data must 
be specified at a given time in order that the equations of motion determine uniquely the evolution of the system." 
(see iQ, p. 147). In formulating NGT as a Hamiltonian system, Clayton (see Q) laid most of the groundwork for 
the study of its initial- value problem. As with most Lagrangian field theories, the field equations of NGT separate 
into a set of initial-value constraints that impose diffeomorphism invariance, and a set of evolution equations. In this 
section, we are concerned with the initial-value constraints of NGT, also known as the Hamiltonian and momentum 
constraints. 

It was demonstrated in the previous section that the momentum constraint can be explicitly satisfied, as in (|^), 
and we need not be concerned with it any further. Meanwhile, the Hamiltonian constraint is a second-order, ordinary 
differential equation for </>, the conformal factor. Unlike the momentum constraint, the solvability of the Hamiltonian 
constraint is by no means guaranteed. Indeed, we can gain some insight by recognizing, as did Wheeler (sec [^ , ^ ) 
the similarity between the Hamiltonian constraint (for a moment of time-symmetry) and a scattering problem in non- 
relativistic Schrodinger mechanics, with the curvature scalar playing the role of the potential. If R^^^^^f is positive 
everywhere, then the Hamiltonian will always have a solution, corresponding to scattering off a potential barrier. If the 
curvature scalar is negative but "shallow" enough, then a scattering solution will also exist. Only if /j'^^f^) jg negative 
and "deep" enough will a bound state begin to form, and hence no scattering solution will exist. This is the type 
of behaviour that will be manifest itself below: there are actually situations where the antisymmetric fields could be 
considered perturbations on a GR background, in the sense that they are by no means large, in which the Hamiltonian 
constraint has no asymptotically fiat solution. Of course, the reason for this behaviour is that the derivatives of 
the antisymmetric variables enter the curvature scalar, along with the variables themselves. Thus, one can choose 
values for the antisymmetric variables that are small in overall magnitude, but that are sharply peaked, causing a 
corresponding "deepening of iJae well" , which in turn causes the Hamiltonian constraint to become unsolvable for an 
asymptotically flat spacetimcu 

In Appendix C, we derive two results that are consequences of the existence of a solution to the Hamiltonian 
constraint, given the approp riat e bou nda ry conditions. In other words, the solvability of the Hamiltonian constraint 
implies that the inequalities (C3) and ( |C4| ) are satisfied. However, it is important to keep in mind that these inequalities 
hold only for massless NGT, when /i = 0. The mass term in (||) makes it considerably more difficult to generate similar 
inequalities for the case of massive NGT. Furthermore, we emphasize that these inequalities are consequences of the 
solvability of the Hamiltonian constraint: to demonstrate that these inequalities actually imply that a solution to the 
Hamiltonian constraint exists is a somewhat more involved process, and will be considered in the future. 

Given an initial data set, {7^^, '(/', j^a}, we can use ( p3[ ) to place a limit on the field variables that will allow us to 
solve the Hamiltonian constraint for an asymptotically fiat spacetime. The procedure is actually one of elimination: if 
we find values of the field variables for which (C3) is not verified, then we may conclude that the Hamiltonian constraint 
cannot be solved for such an initial data set, as the proof in Appendix C concludes that if the Hamiltonian constraint 
is solvable, then (C3) is verified.cl This process is continued, making the initial data set evermore conservative, for 
instance by moving closer and closer to a Minkowski space initial data set, until (C3) is verified. At this point, we 
cannot conclude that a solution of the Hamiltonian constraint exists: rather, we conclude that a solution certainly 
does not exist for less conservative initial data. 

If this procedure is to help us in placing limits on the solvability of the Hamiltonian constraint, then we require 
a (non-trivial) function / that makes the left-hand side of (03) a minimum. This function is well-known [ p8[ : 
= K/ (S^ + r^)]^/^. Here ^ is a constant, and the factor of in the numerator is included to make the value of 
the left-hand side of ( |C3| ) independent of ^; ( |C^ ) is obviously invariant under such a scaling of /. The value of ^ is 
then used to make the value of the right-hand side of ( |C3| ) as large as possible. Note in particular that the maximum 
of / occurs at ^ = ±r. 



As noted in Appendix C, the two inequalities (|C3|) and (C4) are used to different ends. We mentioned above 
that (C3) is used to find a most conservative initi al d ata set {7^^, J^a} for which the Hamiltonian constraint is 



guaranteed to be unsolvable. On the other hand, ( p4D is used to find a least conservative ini tial da ta s et for which 
the Hamiltonian constraint is guaranteed to be solvable. In a sense, the two inequalities, ( |C3| ) and ( |C^ ) are used to 
"sandwich" the limiting initial data set where the Hamiltonian constraint is just barely solvable. 



^In such a case, a solution would most likely exist for a non-asymptotically flat manifold or a closed manifold (as in [Q); we 
are not considering that option here. 

*It is the converse which is not necessarily true: ( |C3| ) might very well be verified, yet the Hamiltonian constraint might not 
possess a solution. 
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In this work, we are considering initial data sets in which 7^^ ~ 1, = 0, and where ip has the form of a pulse: 

V'(r) = ^[r/ro]2e-('-'-°)'/"'. (7) 

The dimensionless parameter A gives the overall amplitude of the pulse, while rg and a are lengths that fix the position 
and the width, respectively, of the pulse; the location of the maximum of the pulse is rpeak = f [''0 + ['"0 + 4(t^]^/^]. 
This particular function was chosen because the Gaussian factor causes both and its derivatives to vanish quite 
rapidly away from rg, while the factor of [r/rg]^ ensures that the Ricci curvature scalar in (|^) remains finite as r — > 0. 
In a sense, the pulse is "semi- localized" about r ~ r-poak, allowing us to study the properties of ngt initial data 
without having to worry about this will affect the boundary conditions. For such an initial data set, the curvature 
scalar is more or less peaked about r — vq. We therefore choose f = rpcak, in order to maximize the value of the 
integrand, pR^^^^l The remaining three free parameters. A, a, and rg, are adjusted in the manner discussed above 



until the inequality in (C3) is no longer verified. For instance, since taking A — reduces the initial data set to that 
of Minkowski space, then we expect to find some maximum value A = A^ax which renders the Hamiltonian constraint 
unsolvable. Of course, this is intuitively obvious: if there is to be a data set for which the Hamiltonian constraint 
cannot be solved for an asymptotically flat spacetime, then we expect that data set to correspond to a strong-field 
region. However, as was mentioned above, it is not only the overall scale of the antisymmetric variables that affects 
the solvability of the Hamiltonian constraint, but also the scale of their derivatives. In particular, smaller values of 
the width parameter a will cause ip to become correspondingly more peaked, and the ensuing increase in its derivative 
could feasibly make the Hamiltonian constraint unsolvable. 

As we approach a region of unsolvability in the parameter space, the mass of the system, as measured by an observer 
in the asymptotically flat region, begins to increase without bound. This is to be expected: a large value of the mass 
would arise out of a spacetime of such great curvature that it can no longer be represented by an asymptotically 
flat manifold. It is therefore of great interest to measure the mass of the system as a function of the parameters of 
the pulse defining ■0- Through an analysis of the surface terms that arise out of the canonical decomposition of the 
action , it can be shown that the ADM mass (see ) for asymptotically flat initial dataCl is 

M = M-^i> 7"''Va[(/.] d56 (8) 

where 9E is the boundary of the spatial slice and where M is the ADM mass of the conformal metric and M is the 
ADM mass of the physical metric. For the conformally flat field variables we consider, and given the fall-off that we 
have assumed, it is a simple matter to demonstrate that (|^) reduces to 

M = -2 I l^^[(t)\r'^dr (9) 



where the derivative operator A^[] is defined in (CI). Once a solution to the Hamiltonian constraint is found, it 
can be inserted into (^), and the mass M calculated; note that this mass is guaranteed to be positive as and 
the Bonnor term are negative definite. Alternately, if we assume that (j) ^ 1 + M/2r asymptotically, where M is a 
constant, we can read off the value of M by inspecting the behaviour of (p in the asymptotic region. 

In the next section, we discuss a numerical scheme that can be used to analyze the initial-value problem of ngt. 
In addition to presenting a typical solution to the Hamiltonian constraint, we use this scheme and the conditions 
discussed above to place limits on the solvability of the initial-value problem in not. 

IV. NUMERICAL ANALYSIS OF THE INITIAL- VALUE PROBLEM IN NGT 

As written, the Hamiltonian constraint ( ph| ) is a second-order, ordinary differential equation for the conformal factor 
(j). There are two important points to note. Firstly, suppose we choose GR initial conditions: i.e., 7^^ ~ 1, R = r 
(standard spherical coordinates), and vanishing antisymmetric-sector data (■0 = and = 0). The Hamiltonian 
constraint then reduces to r^9^[0] -I- 2rdr\(j)\ = dr[r^ dr\(t)\\ = 0, whose solution is = 1 + M/2r, where M is a constant 
of integration, physically interpreted as being the mass defined by (m. Secondly, when /i 7^ 0, the Hamiltonian 



^For our purposes, asymptotically flat initial data is defined by the usual conditions on the symmetric-sector functions (as in 
GR), and "faster" fall-off in the antisymmetric-sector variables. This is defined explicitly in fesl, p. 116. 
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constraint is a non-linear, second-order, ordinary differential equation for the conformal factor, in contrast to GR, 
where the Hamiltonian constraint for a vacuum spacetime is linear in the conformal factor when the extrinsic curvature 
vanishes (moment of time symmetry); to see the non- linear behaviour, matter or energy of some kind must be added 
to the system. Of course, this non-linearity in the Hamiltonian constraint is a manifestation of the fact that, loosely 
speaking, the antisymmetric-sector fields behave like matter or energy to the symmetric-sector fields, albeit with a 
nontrivial interaction. Fortunately, since the Hamiltonian constraint in spherical symmetry is an ordinary differential 
equation, we can use a fairly simple scheme for its solution, and it is this procedure that we describe in this section. 

Since we are considering initial data that are concentrated in an area of the initial slice, in those areas of the 
initial slice where the initial data of the antisymmetric sector are vanishing or nearly so, will be more or less 
approximated by its GR solution, cj) ^ 1 + M/2r, where M is either some constant or a linear function of r. The first 
case corresponds to a Schwarzschild-like region of the initial slice, while the second case is physically equivalent to a 
flat, Minkowski region of the initial slicefl We would expect to find the former behaviour in the asymptotic region, 
while the latter behaviour would be prominent near the origin, given our choice of tp. The only region which remains 
in question is the intermediate region, where the antisymmetric-sector initial data are non-trivial; however, this 
essentially becomes a region of transition that serves the purpose of matching the Minkowski-type solution near the 
origin to the Schwarzschild-type solution of the exterior. The obvious generalization is therefore to assume that AI no 
longer has this simple behaviour, and to make the substitution M M{r). We will therefore use (f){r) — 1 + A/(r)/2r 
for all r, and solve the Hamiltonian constraint for the function M(r). 

Upon substituting tp — \ ^ M(r)/2r into the Hamiltonian constraint (pb|), we obtain 



2r 



rdr[M]) + \ 



1 f cosi/' 
^11 



"(j^s)^ fJ^'^rcj}^ sin^ tp 



27110^ 



I67" 



(10) 



where (j) is to be considered shorthand for (j) — I + M{r)/2r. This equation is finite-differenced across the spatial grid, 
using centred differences for both the first- and second-order derivatives. 



dr[M] 



M,.- 



2Ar 



and 



d^M] - — 
or,; 



Sri 



M,+i - 2M, + M,_i 
Ar 



where Ar is the (uniform) spatial grid spacing. El This yields an equation of the form A/^+i — 2Mi + Mi-i = f{M), 
where the notation f{M) is used to represent the fact that the right-hand side has some functional dependence on 
the value of M at grid, points i + i, and i — 1, along with the other variables. An implicit relaxation scheme is used 
to solve this equationD: beginning with an initial guess for the Mi, we write M^^-^ — 2M| -I- M^_-^ = f{M^^^) and 
solve for the This process is repeated until a precision criterion is reached: in our case, we chose to calculate 
the integral of the left and right sides of (|l^) for a given iteration. If the difference between these two is smaller than 
some predetermined value, then the solution has converged and the iteration ceases. The accuracy of the results will 
be discussed at the end of this section. 

Evidently, ( [To| ) and its finite-differenced counterpart are to be augmented by boundary conditions on M that will 
fix the exact physical situation that is being studied. In this work, we have chosen to take M to be linear in r near 
the origin. Numerically, this boundary condition is most easily implemented by noting that if M is linear in r, then 
the average value of M at a given point coincides with the value of M itself at that point: {fJ,M)i = Mi. Evaluating 
this at the origin, where M = 0, we have (^M)o = 0, where we have labeled the origin by gridpoint i = Q. In the 
asymptotic region, we have taken M to behave like a constant, and thus its derivative must vanish: cJriA-/] = 0. This 
guarantees that — > 1 as r — s- +00. Numerically, this is encoded using either a backward difference based on the grid 
point n, or a forward difference based on the grid point n — 1, where n labels the outer-most grid point: because the 
derivative is made to vanish, both formulations are equivalent. Therefore, in the asymptotic region, we take either 
(vM)„ = Oor (AM)„_i =0. 

In Figure |l|, we give a plot of the conformal factor of massless not for the initial conditions 7^^ = 1 and = 0, 
with tp given by (R) with A — 0.7 and rg = 40 in units where a = 10; in this and in every other case in this paper, 
the units of r, rg, cr, and M will be assumed to be the same. The solution was found using the implicit relaxation 



^If M is a linear function of r, then M/2r is some constant, and therefore is a constant other than unity. However, by 
rescaling the basis vectors, we can absorb this constant and make <j) = 1, just as in the case of Minkowski space. 

^The notation is standard: see, for instance, p. 48, for definitions of the centred-differencing operator, 5, and the averaging 
operator, as well as the forward- and backward-differencing operators, A and V, which will be used later. 

®Such schemes are described in the standard references; see, for instance, |4ll], p. 86. 
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FIG. 1. The conformal factor of massless ngt as a function of radial distance from the origin, for the initial conditions 
7^^^ = 1, tp = A[r/ro]^e~'' ^'^ , and j^^ = 0. We have taken A = 0.7 and ro = 40 in units where a = 10. The radial variable 
r is measured in the same units as ro. These data have an adm mass of M « 22.340 in the same units as r. 




FIG. 2. For the same initial data and parameters as Figure |l|, we overlay a magnified region of that graph (lower curve) onto 
a plot of the GR conformal factor, 1 + M/2r, where M ~ 22.340 in the same units as r. 
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scheme described above. The values of the parameters were chosen to place the simulation in an intermediate region 
of the parameter space: sufficiently far away from Minkowski space so as to avoid making the simulation trivial, and 



sufficiently far away from the regions of unsolvability discussed in §111 to avoid infringing on some limiting case that 
could raise questions about the stability of the numerical solution. Throughout, we have chosen rmax, the location 
of the outermost grid point, to be at least Tpoak + lOcr: this guarantees that '!/'('"niax) ~ and that we can impose 
asymptotically flat boundary conditions. The grid spacing was taken to be Ar = 0.015 in the same units as r. In 
other words, the solution depicted in Figure |l| represents a generic conformal factor for the initial data set used in this 
work, with the antisymmetric sector variables being strong, but not pathologically so. The ADM mass of the system 
(see the discussion in §111, in particular (||) and (^) for the definition of the ADM mass in NGt) is calculated to be 
M « 22.340. The graph of Figure || overlays the GR conformal factor 1 + M/2r, where M ft! 22.340, atop a magnified 
section of the conformal factor generated from the numerical code. This graph demonstrates that M{r) goes to a 
constant in the asymptotic region, despite the large width (cr — 10) of the initial data. 

The solution displayed in Figure |l| evidently possesses many of the features described earlier. In fact, the radial axis 
is clearly divided into three regions, corresponding to the central peak of ■0, and the regions a few widths a above and 
below this peak. For our choice of parameters, the central peak of ip is located at rpcak ~ 42.36. The region < r < 30 
corresponds to a Minkowski space solution, while the region r > 60 corresponds to a Schwarzschild-type solution, 
= 1 + M/2r, where M is a constant; in this case, M « 22.34. In the latter region, which is at least two widths 
away from the central peak at rpcak, the value of "0 is, by construction, entirely negligible. The antisymmetric sector 
field variables behave essentially like an shell of energy to the symmetric sector field variables, which then collapse to 
their GR form: i.e., Schwarzschild-like. This much is hardly surprising, and indeed, could easily have been predicted 
based on a perturbative analysis of the Hamiltonian constraint. However, what is somewhat more surprising is how 
quickly the solution to the Hamiltonian constraint relaxes to its GR form outside the influence of the antisymmetric 
sector field variables. In fact, it is clear from Figure ^ that in the region 50 < r < 60, has essentially returned to 
its GR form: the difference between the two curves in the graph is at its greatest less than 0.5%, and this despite the 
fact that this region lies within two widths of the central peak of V'- Indeed, at r = 50, '0(5O)/V'(rpcak) ~ 0.54, which 
is hardly negligible. This demonstrates the rather surprising result that the antisymmetric field variables can very 
well be non-trivial, and yet the gravitational system as a whole still behaves very much like GR, at least as far as the 
initial- value problem is concerned. 

In Figures | through || we present plots of the measured mass of the gravitational field versus the various parameters 
of the pulse; the value of the parameters that remain fixed for the plot are given in the figure caption. For the massless 
NGT calculations, we also supply in the figure caption the bounds on the parameters determined from ( |C^ ) and (|C4|); 
note that our initial data allows us to use the flat-space Sobolev constant. Of particular importance is the behaviour 
of the mass as a function of the pulse width, cr, for the case of massive NGT. We see from Figure ^ that there is a 
maximum value of a above which the Hamiltonian constraint is no longer solvable. Loosely speaking, this implies 
that the pulse (^) must be sufficiently small and sufficiently peaked (for a given amplitude) in order for a solution to 
the Hamiltonian constraint to exist. This seems somewhat counter-intuitive; however, when /i 7^ 0, we see from (|^) 
that the term proportional to contains a factor of sin^ tp, and is also proportional to 0^. An increasing value of 
the width causes the factor of sin^ -0 to become non-trivial in a larger region of the spatial slice, thereby increasing 
the energy contained in the gravitational system. Consequently, this forces to stray from its assumed — s- 1 fall-off, 
and the Hamiltonian constraint becomes unsolvable. It therefore stands to reason that this highly non-linear term 
can affect the resulting solution to a great extent. Of course, a similar behaviour can be observed in the abscence of 
this term: loosely speaking, increasing the width increases the amount of antisymmetric field that is present on the 
spatial slice, therefore adding to the energy of the system. However, we have found that this occurs at width values 
on the order of cr « 217, where M « 3190. The inequalities in Appendix C give 51000 ^ cr^ ^ 55000. We conclude 
that the non-linear term in (|^) simply compounds this problem. 

We note from the figures that there are clearly defined asymptotes in the mass as a function of the parameter being 
investigated; these asymptotes correspond to the divergence of the mass as we approach a region of unsolvability of the 
Hamiltonian constraint in the parameter space. In particular, the mass plots for massless NGT have asymptotes which 
respect the bounds set by the results of Appendix C. For instance, from Figure || we see that the asymptote of the 
mass as a function of the amplitude lies at A « 0.70. This compares well with the bounds oi A ^ 0.51 and A ^ 0.81 
set by (p4) and (C3), respectively. Similarly, the asymptote of the mass as a function of the width squared (see 
Figure ^)hes at about cr-^ » 0.005 for tq = 40. This is to be compared with the bounds of cr^ ^ 0.0046 and a'^ ^ 0.090 
obtained from ( |C3| ) and (C4), respectively. Finally, the asymptote of the mass as a function of the position of the 
pulse rp ( see Figure 0) lies at tq » 620 for cr^ = 1.6, in comparison to the values of tq ^ 160 and tq ^ 750 obtained 
from (C4) and (|C3|), respectively. 

Generically, the asymptotes of the mass are more clearly defined for massless NGT than for massive NGT. This is 
most pronounced when plotting the mass M versus the position of the pulse tq : it is clear from Figure |^ that the 
asymptotic value of tq has not been reached. The reason for this is the same as mentioned above: the non-linear 
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FIG. 3. The mass of the gravitational field versus the amplitude A of the pulse for massless ngt, where cr^ = 3 in units 
where ro = 40. The inequalities in Appendix C give 0.51 ^ ^ ^ 0.81. 



M 




FIG. 4. The mass of the gravitational field versus the amplitude A of the pulse for massive ngt, where cr^ = 3 in units where 
ro = 40. 
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FIG. 5. The mass of the gravitational field versus the squared width of the pulse for massless ngt, where A = 0.15, and 
where ro = 40 in the same units as a. The inequalities in Appendix C give 0.0046 ^ cr^ ^ 0.090 for the left asymptote, in those 
same units; the leftmost point on the graph is at cr^ fa 0.0060, where M w 566. 
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FIG. 6. The mass of the gravitational field versus the squared width of the pulse for massive ngt, where A = 0.15, and 
where ro = 40 in the same units as a. 
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FIG. 7. The mass of the gravitational field versus the position ro of the pulse for massless ngt, where A = 0.15, and where 
cr^ = 1.6 in the same units as ro. The inequalities in Appendix C give 160 ^ ro < 750 in those same units. 
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FIG. 9. The mass of the gravitational field of massles NGT as a function of the outer grid point, rmax, for the initial conditions 
7" = 1, V = A[r/ro\^e~''''~''°^ , and fs = 0, where A = 0.7, and with a = 10 in units where ro = 40. 

TABLE I. The residual Hamiltonian He and the ADM 
mass M as a function of the grid spacing Ar for mass- 
less NGT, for the initial data 7^^ = 1 and j^s = 0, with 
^ = A[r/ro]^e-('-''o'''/'"', where A = 0.7, and ro = 40 in 
units where ct = 10. 



Ar M He 

0300 22.317 4.525x10"^ 

0.150 22.334 1.129x10"^ 

0.075 22.338 2.815x10"^ 

0.030 22.340 4.554x10"'' 

0.015 22.340 1.356x10"'' 



term in (||) begins to grow significantly, and the numerical code fails to find a solution. It is also inte resting to note 
that when determining the bounds on the region of solvability for variations in a (see Figure I), (jcl) seems to give 
particularly bad results, being ofi^ by a full order of magnitude. Although the predicted bound does not contradict the 
actual result, and no claim was made as to the accuracy of (|C^), the reason for why there is such a large discrepancy 
is unknown at this time. 

When calculating the ADM mass of a system, we must take into account the fact that our numerical grid is not of 
infinite size, whereas (||) implicitly assumes that the evaluation is performed at infinity. In Figure ^ we show how the 
ADM mass varies as the size of the grid is increased, for the initial data = 1 and j^s = 0, with ip given by (|^), 
where A — 0.7, ro = 40 and a = 10. As expected, the mass quickly converges to a definite value once the grid is 
sufficiently large. In our case, the mass has converged to within 0.01% when r^ax — 70. At this point, ip has dropped 
to approximately 0.04% of its peak value. This justifies our choice of rmax = fpcak + lOcr. 

Finally, it is important to verify that the solution (f) obtained by numerical means actually satisfies the Hamiltonian 
constraint, and if so, how well. This is somewhat of a thorny issue since, for a vacuum spacetime, the correct solution 
must make the Hamiltonian constraint vanish. It is therefore inheritly difficult to judge how well the solution satisfies 
this criterion without imposing some arbitrary precision level. One possibility, similar to that discussed by Bernstein 
et. al. (see Appendix A. 3 in ||4^) is to use the fact that in a non- vacuum spacetime, the Hamiltonian must be 
equated to the matter contributions: 7i = 2pp. Thus, the error in the determination of <j) can be viewed as a "matter 
contribution" to the right-hand side of the Hamiltonian constraint. Consequently, the mass of this contribution can be 
calculated and then compared with the known mass of the hypersurface, which then gives a measure of the precision 
of the solution. 

In order to verify that the solution satisfies the Hamiltonian constraint to a sufficient degree, we calculate the 
integral of absolute value of the Hamiltonian constraint: 
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He^ J \n\d^x, 



where the integration is performed over the volume of the shce S. We refer to i?e as the residual Hamiltonian, as it is 
a measure of that part of numericaUy-determined solution that violates the constraint. The residual Hamiltonian is 
then compared with the ADM mass M of the hypersurface. This process is repeated for successively smaller values of 
the grid spacing, until the ratio H^/M is reduced to the desired level of precision. In Table | we list the results of such 
a calculation in the case of massless ngt for the initial data 7^^ = 1 and = 0, with -if) given by (Q) where A = 0.7, 
rg = 40, and a = 10. The ADM mass and the residual Hamiltonian are given in the same units as Ar, the grid spacing. 
We see that as the grid spacing is reduced, the ADM mass converges to a definite value, in this case M w 22.340, 
and residual Hamiltonian decreases roughly as (Ar)^. This is in keeping with the fact that our discretization of the 
Hamiltonian constraint was valid to second-order in the grid spacing. 



V. CONCLUSIONS 



We have described herein the foundation necessary for performing a thorough numerical study of spherically sym- 
metric systems in NGT. By virtue of its construction, this formalism is not limited to a static gravitational field; this 
is particularly relevant for the spherically symmetric, Wyman sector initial data we have described in this work. In 
fact, the antisymmetric field variables are true dynamical degrees of freedom, not limited to being diffeomorphically 
related to a physically-equivalent static system by some analogue of Birkhoff 's theorem. 

We have demonstrated that the initial- value problem of NOT can be solved by numerical means, and that the 
solutions behave essentially as one would expect, based on a careful physical analysis of the situation, as described at 
the beginning of §[V. However, for a given initial data set, we cannot always solve the initial- value problem. We were 
able to place limits on the solvability of the Hamiltonian constraint for a particular distribution of the antisymmetric 
field variables. In itself, it is not particularly surprising that the initial- value problem cannot always be satisfied for all 
possible field configurations. What is surprising, however, is that certain seemingly acceptable configurations do not 
form a well-posed initial- value problem. In particular, for the case of the massive theory (/.t ^ 0), it was found that 
the antisymmetric field variables had to have particularly rapid fall-off in the asymptotic region in order to guarantee 
that the initial slice be asymptotically flat; otherwise, the final term in (pb|), proportional to both svtiip and <^^, would 
feed back into the solution and destroy the desired asymptotic behaviour. 

The behaviour of the ADM mass, or rather its not equivalent, was found to be a good indication of the solvability 
of the Hamiltonian constraint. Indeed, far from regions of unsolvability, the ADM mass was very small, corresponding 
to a weak field regime. On the other hand, near regions of unsolvability, the ADM mass became very large. But in the 
latter cases, the mass always began to display signs of reaching an asymptote. Again, this corresponds well with the 
underlying physical picture: a large mass cannot be obtained from an asymptotically flat spacetime, so demanding 
that the mass be large by increasing the size of the antisymmetric field variables causes us to step off the constraint 
surface, thereby violating the Hamiltonian constraint. 

Of course, there remains much to be explored. For instance, we have not touched upon the topic of the Cauchy, 
or evolution problem in NGT. Once we are satisfied that a given field configuration properly satisfies the initial-value 
problem, we can then use the NGT evolution equations to determine the evolution of the system. The definition and 
existence of apparent horizons are also pivotal to the structure of NGT; this will be dealt with in a future publication. 

Beyond this, we could modify the formalism (and hence the numerical code) to accommodate matter sources, 
thereby moving beyond the purely vacuum case that we have been considering here. Ideally, matter fields would be 
combined with an evolutionary approach, and we could make some headway into determining precisely what it is that 
happens in NGT when a stellar object expends its fuel and begins to collapse. 

Perhaps the most critical aspect of NGT is that it remains as a candidate for a non-singular theory of gravitation. 
These claims are based on the study of a single, spherically symmetric, vacuum solution to the field equations of 
massless NGT (see and |^^) as well as a certain amount of indirect theoretical evidence (see |44 ). Meanwhile, 
some workers have claimed that NGT does exhibit the same singularity properties displayed in GR based on perturbative 
analyses of the field equations (see |4|]). Ideally, we would prefer to consider initial data that describes a stellar model, 
essentially along the lines of the classic Oppenheimer-Snyder collapse problem in GR (see [Q), and "let the thing 
go" . However, the mathematical difficulties involved in performing such feats has for some time relegated the topic 
to the area of pure speculation. Given the formalism we have outlined herein, this ideal is much more attainable than 
it once was. 
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APPENDIX A: THE SPHERICALLY SYMMETRIC FIELD EQUATIONS OF NGT IN FIRST-ORDER 

FORM 



A detailed description of the Hamiltonian form of NGT can be found in | |2^ . In this appendix, we give a brief overview 
of the reduction of these field equations for spherically symmetric systems. In particular, we will be considering 
exclusively the "Wyman sector" of NGT, also known as the electric sector (see 19 2^,^). The inverse fundamental 
tensor will therefore take the form g^^ = e±^e± — 7°''ea ^ eb, where the basis vectors are defined by Ne± = dt — N°-da 
and ea = da 

The spherically symmetric form of a nonsymmetric tensor quantity can be found by a Killing vector analysis 
(see [p3|, p. 94). The spatial part of the inverse fundamental tensor and the extrinsic curvature are 



.ab 



1 





^22 7(23]/ sin ^ 
-7123] /sine 722/ sin^^ 



and 



Kab 










-j[23] sin^ 





J[23] sin^ 
K22 sin^ 



notice how we reserve the indexed symbols for the tensorial quantity with factors of s'mO removed. The volume 
element is written 



p sm ( 



smf 



The lapse function depends on the radial variable (and the time), TV = N{r,t), while the shift vector reduces to a 
single non-trivial component: |A^°| [A^^,0,0]. The vectors jp"! and \ Wa\ vanish by definition, being both native to 
the "Papapetrou" or magnetic sector of NGT. The non-zero components of the fundamental tensor are then 



Gn = (7")-\ 
We define the following: 



7 



22 



(722)2 ^ (7(23] )2' 



and G 



[23] 



7 



[23] 



(722)2 _,_ (7[23])2 



722^^22+ 7''^' J[23], and j23=7[23]i^2 



22 ■ 
7 J[23]- 



For spherically symmetric Wyman sector initial data, the non-zero Lagrange multipliers and connection coefficients 
can be solved for, to wit: ai = dr[hiN], and — 7^^ai, while 



Tl, = i9.[ln(p2[(7 



22n2 



ia„[ln(p27n[(722)2 + (^[23])2]l/2)]^ 



and 



^22 — -7"^MG22r22 - G[23]Ai3]. 

The non-trivial components of the tensor |A^^| are 

A?3 = ^G227"'l(a.[ln(p7'''l) - a,[ln(p722)]), and A^g = 7"(G[23]r?2 + G22A23 



Finally, the components of the NGT surface Ricci tensor are found using the connection coefficients and the components 
of |A|J^| given above: 



nNS(3) _ 
^11 — 



-25.[r22] + 2r22(r}i-r22)-2(A23)2 



R 



NS(3) _ 



"9i-[r22] + r22rii + 2A23A23 -1- 1, 



and 



NS(3) 
[23] 



9r[A23] 



22^13- 
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APPENDIX B: THE CONFORMALLY-TRANSFORMED FIELD EQUATIONS AND CONSTRAINTS 

FOR SPHERICALLY SYMMETRIC NGT 



In order to approach the initial- value problem in NGT, in particular so as to study the solvability of the Hamiltonian 
constraint, we take the same route as in GR, which consists of conformally scaling the field variables by a factor 0, 
thereby transforming the Hamiltonian constraint into an equation for this conformal factor (see ^^). The variables 
|7ab| and \Kab\ are chosen to transform in some particular way, which then determines the transformation rules for 
the connection coefficients. Although this can be done in general in GR, the process is tantamount to inverting the 
compatibility condition, which in NGT is not nearly as trivial as in GR (see, for instance, p8| , ^ ). 

We choose |7ah| and \Kab\ to transform as ^ab — > 4>^lab and Kab — > 4' '^Kab (see Q); from this we determine that 
ryab _^ (j)-i^ab _^ (p^Gab, while K'^b (ji-^K^b and ^ (p^^fb- Note that we choose the transformation 

rule for the antisymmetric components of the fundamental tensor so that the determinant of the fundamental tensor 
transforms homogeneously m (j): p ^ The non-trivial Lagrange multipliers are found to transform as — > (f)~'^a^ 
and ai ai, while the connection coefficients transform as 

Tn ^ -f29^[ln0], Vj^ ^ Vj^ + 2dr[\n(j)], and ^ - 2-f^'^G22dr[ln(j)]. 

Finally, the components of the tensor |A^^| transform as 

A?3 ^ A?3 and A^g ^ A^g -I- 27"G239r[ln</']. 

Despite its complexity, the transformation rule for the momentum constraint is relatively straightforward: 

Hi^Hi- 8p[K\ + 2if22]a,[ln0] - -[K\ + 2K\]dr[p] - i^nS^H"] + 2K22dr[p7^^] + 2j[23]dr[pi^^^^ 

+ W^dr[K22] + 4p7'''l5,[j[23]] - Sp[K\ + 2K\]dr[lii4>] = 0. 

Meanwhile, the various terms in the Hamiltonian constraint transform differently, yielding 

n = sp^'^m^] - 8p</-a.[0][7"rli + 2^''rl, + 271^^1 Af^a,] - p^P'r""'^'^ 

= 8p^A,[0] - p02i?^s(3) _ 2p4>-MK\f + 2K\K\ - - \^?pct>'l^''^G^23} = 0, 

where in the last line, we have used the definition ( |Cl| ) to write the Hamiltonian in the familiar Lichnerowicz form. 

The evolution equations for the components of the fundamental tensor can be shown to transform in a somewhat 
simple manner: 

9J7"] = iV^^.b"] - 27"a,[7Vi] - 47VS"a,[ln0] ~ 2N cj^-^^^^f K^^, 
5, [7^2] = dr[NW^] - ANW^dr[\ncP] - 27V0-6 [(7^2)2 + [^^^^\f]K22. 

and 

a47[23l] = a.iiVS''''] - 47V^7'''l5.[ln<^] - 2Nr^[{^^^f + 

The evolution equations for the components of the extrinsic curvature would be equally as simple, if not for the 
complicated transformation rule of the components of Zab] in fact, we find that 

dt[K^^] = N^dAK^x\ -h 2K^^dr{N^\ - 27Vii^n5,[ln0] - TV^^Zn - \[?f^<t>''Grr^'^''''^G^23\. 

dt[K22] = N^dr[K22] - 2N^K22dr[ln<j)] ~ N<j,^Z22 + iAi'Af</''G227''^'G[23], 



and 
where 



and 



5t[j[23]] = ^'a,[j[23]] - 2N'j[23]dr[\n<j>] - Ncl>^Z[23] + N cj,'' [-f^^^^G 22 f " G[23]] 

Zn = -4[(/.92[0] _ (9,[0])2] + 2[a,[ln7V] + 2rii - 2Tl^]cl)dr[<l>] 

+ (^'^ + a,[ln7V](r}i - a,[ln7V]) - ^^[InTV]]^^ + [(X\ + 2K\) - 2K\\K^^ct^-\ 
Z22 = -27"G22[R'[0] + (drW\f\ - 2[G22a,b7"] - +P7"G229r[lniV]](/.a,[0] 

+ + dr^^N\V\^\<\? + \K\K22 + 2^23/3]'/'"', 

Z[23] = 27"G[23][09,'[0] + {drW\f\ + [2G[23]5.[P7"] + 2A^23] + 2p7"G[23]a,[ln7V]]0a,[(/.] 
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APPENDIX C: TWO CONSEQUENCES OF THE SOLVABILITY OF THE HAMILTONIAN 
CONSTRAINT FOR MASSLESS NGT IN SPHERICAL SYMMETRY FOR A MOMENT OF 

TIME-SYMMETRY 



The Hamiltonian constraint of massless NGT for a moment of time-symmetry can be written in the form 8(/)A^ [0] 
^2^NS(3)^ where the derivative operator Aj[] is defined by 



(Cl) 



We derive herein two criteria that are consequences of the existence of a function (jj, solution of the Hamiltonian 
constraint of massless NGT, having the appropmate properties (0 > and — 1 —> as in the asymptotic region 
in order to guarantee finiteness of the energy. )□ 

Let / be a twice differentiable function with asymptotic fall-off. Since > 0, then / can be written as the product 
f = (f)u, where u is also twice differentiable and has the same fall-off as /. Hence, 



(a^[/])2 ^ (t)^{dr[u])^ + 2(f>udr[(f>]dr[u] + {dr[(t)])^ = <l>^ {dr[u]f + 09^ [(/)] 9^ [u^] + {dr[(t)])^ . 

Multiplying both sides by P7^^ (to form a scalar density) and integrating gives 



(C2) 



psm9j''{dr[f]rd''x^ / sin6l(p7^V - V5r[P7 ]5,.[</'])d''a;= / psin6l(7^V (^^M) - 0m'A^[<?!)]) d^x, 



where the cross term ipdrlcjildrlu^] in ( |C2| ) has been integrated by parts; the resulting surface terms vanish due to the 
fall-off of u and (j). Now, the first term on the right hand side is manifestly positive, so that we may conclude that 



psmeY\dr[f]y d-^x > - J psmeu''(j)A^[(j)]d-^x. 

However, assuming that (/> is a solution of the Hamiltonian constraint and using the fact that (j)u — /, we may write 

1 



psin9-f''{dr[f]yd-'x > - 



psm 



(C3) 



This condition must be satisfied for all functions / with asymptotic fall-off. Alternately, if the Hamiltonian constraint 
admits a solution cj), then the NGT scalar curvature is such that all functions / having asymptotic fall-off will 

satisfy Q. 

For any asymptotically flat Ricmannian three-manifold, the Sobolev inequality ]5l| ] states that there exists a positive 
constant C such thatEj 



psin9Y^{dr[C]) d'^x ^ C / psmOC d'^x 



1/3 



for any infinitely differentiable function C having compact support: ( G . For a flat manifold, the Sobolev constant 
is C = j2'^/^tt'^/^ . For a non-flat manifold, the value of C is most probably smaller [Q. Using the Sobolev inequality 
and the Holder inequality (see p. 52), 



psineCR^^^^'' d^x / psin6'C^ d^x 



1/3 



psm 



^|^NS(3)|3/2^3^- 



2/3 



we can transform (C3) into 



psin0|i?NS(3)|3/2^3^- 



2/3 



^ 8C. 



(C4) 



^This proof is a slightly modified version of that which appears in the notation is similar, although it has been modified 
when necessary. It has been modified to fit your screen. 

^"Note that our constant C diflfers from that of Cantor and Brill; see, for instance, (10) in In fact, where we use C, they 
use C'\ 
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The two inequalities (C3) and (C4) give qualitatively different results. The first of these is used to demonstrate 
that a data set does not support a solution of the Hamiltonian constraint, while the second is used to prove that a 
data set does support a solution of the Hamiltonian constraint. 
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